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Abstract: We study combined leptoquark {IQ) single and pair production at LHC at 
the level of detector simulation. A set of kinematical cuts to maximize significance for 
combined signal events has been worked out. 

It was shown that combination of signatures from single and pair production not only 
significantly increases the LHC reach, but also allows us to give the correct signal interpre- 
tation. In particular, it was found that the LHC has potential to discover IQ with a mass 
up to 1.2 TeV and 1.5 TeV for the case of scalar and vector iiQ, respectively, and IQ single 
production contributes about 30-50% to the total signal rate ioi IQ — I — q coupling, taken 
equal to the electromagnetic coupling. 

This work is based on an implementation of the most general form of scalar and vector IQ 
interactions with quarks and gluons into CalcHEP/CompHEP packages. This implemen- 
tation, which authors made publicly available, was one the most important aspects of the 
study. 
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1. Introduction 



Boson fields mediating lepton-quark interactions naturally appear in various extensions of 
the Standard Model which is known to be theoretically incomplete. Leptoquarks appear 
in the framework of Grand Unified theories (GUT) where quarks and leptons are unified 
in one matter multiplet [|^, in the SUSY models with R-parity violation (in this case the 
mediator of the lepton-quark interaction is a squark or a slepton) , as well as in composite 
models of leptons and quarks . 

In general, boson fields give rise to violation of the baryon and lepton numbers (lead- 
ing to fast proton decay) and flavor changing neutral current(FCNC) processes which 
are strongly constrained by experiment. In principle, those bosons should be heavy (~ 
Mpianck ~ 10^^ GeV) to suppress these unwelcome processes. On the other hand, fast 
proton decay and FCNC problems can be avoided if the boson mass even is of the order 
of the electroweak (EW) scale. In the case of leptoquarks {SQ), which exclusively induce 
lepton-quark interactions, the problems above are solved: 1) the fast proton decay problem 
is absent, since leptoquarks conserve the baryon and lepton numbers (in the case of squarks 
mediating lepton-quarks interactions, in R-parity violating SUSY models, the proton decay 
is preserved when only one, lepton or baryon number is conserved); 2) the FCNC prob- 
lem is also absent under the assumption of flavor diagonal form of leptoquark-lepton- quark 
(BQ — I — q) interactions. 

Numerous experimental IQ searches at HERA (e.g. Q) and at the Tevatron (for recent 
results see, e.g. Q and references therein) gave no positive results so far and provided only 
limits on IQ masses and IQ couplings to leptons and quarks. Eventually, the CERN Large 
Hadron Collider (LHC) will be able to extend significantly the reach for IQ masses and 
couplings (see, e.g. [^, |8|, ^J). In [10, 11| it was shown that Next-to-Leading-Order 
(NLO) corrections to IQ pair production are positive and non-negligible ( up to 20% at the 
Tevatron and up to 90% at the LHC) and lead to further extension of the collider reach 
in IQ searches. Leptoquarks can be produced not only in pairs (from gluon splitting) but 
also as a single particle in association with lepton |l^, ^ . In ||l^ authors noticed the 
importance of the combination of IQ single and pair production analysis for the Tevatron 
collider, which partially motivated the present study. 

In this article, we perform a new detailed study of IQ production and decay at LHC. 
There are several motivations for this work. Firstly, we stress the importance of the study 
of combined IQ single and pair production, since both contribute to the same signature 
at the detector simulation level. Therefore, combination of signatures from IQ single and 
pair production not only significantly increases the LHC reach, but also allows us to give 
the correct signal interpretation. Secondly, we study the complete set of IQ interactions 
including scalar and vector IQ production and decay. Finally, one should stress that this 
study is based on the implementation of the complete IQ model including the most general 
form of scalar and vector IQ interactions with quarks and gauge bosons (including gluons) 
into CalcHEP/CompHEP packages 2C]. This is one of the most important aspects of 
this work. Such an implementation allowed us to study the effect of vector IQ interactions 
with gluons via anomalous couplings of the most general form. 
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This paper is organized as follows. In section 2, we describe the general effective 
Lagrangian used in our study, as well as the implementation of this Lagrangian into the 
CalcHEP/CompHEP software package. In section 3, we present signal rates for IQ single 
and pair production at the LHC for the case of scalar and vector E^. In section 4, we 
perform an analysis of signal versus background at the detector level. Finally, in section 5, 
we draw conclusions on LHC potential for leptoquarks search. 

2. ]Q model and its implementation into CalcHEP/CompHEP packages 
2.1 The model 

Following ll^, |l^, we use an effective Lagrangian with the most general dimensionless, 
SU{3) X SU{2) X C/(l) invariant couplings of scalar and vector leptoquarks to leptons and 
quarks with lepton and baryon number conservation: 

^ = 4f|=0 + ^|F|=2 + ^'' (2.1) 

Lagrangian ^^p^-Q2 describes Yukawa type interactions of with leptons and quarks 
(IQ — I — q), changing the fermion number F by or 2, respectively, where F = 3B + L, 
B is the baryon number and L is the lepton number. ^ conserves the baryon and 

lepton numbers and has a flavor diagonal form: 

+ {hiiqi^L + hiRdR'y''eR)Ui^ 

+ hiRURj^'eRUi^ + HLqLf-f^^LU^t, + h.c, (2.2) 



^^\F\=2 = {9iLqliT2h + giRujieR)Si 
+ gmdReRSi + g?,LqliT2T(.LSz 

+ {92LdRl^iL + 92Rqll''eR)V2^ 

+ 52n^7^^LV2/. + h.c. (2.3) 

where Tj are the Pauli matrices, gz, and II are SU(2)l quark and lepton doublets, re- 
spectively and ur, dR, and cr are corresponding singlet fields; charged conjugated fields 
are denoted by = C/^. We follow the classification from [|l7|. Table || summa- 
rizes the complete set of scalar and vector BQ fields appearing in Eq. (^)- (|2.3|) : Si, Si, 
S3, R2, R2, and V2, V2 1 U^, C/f, C/3, respectively. Scalar {Si, Si, S3, R2, R2) and vector 
(Vg^, Vg^, C/f, f/f, C/3) IQs have Yukawa-type couplings to quarks and leptons denoted by 
{giiL,R), giR, gsL, h2(L,R),h2L) and {g2(L,R),g2R,h^L,R)MR,h3L) respectively, in Eq. (^- 



( |2.3| ). Couplings appearing in the Feynman rules of BQ — q — I {BQ — q — v) interactions 
(which we denote by Xiilq), ^Rilq) and Xiivq)) are trivial linear combinations of g, h,g, h 
couplings presented in Table |l|. The table also presents a notation for names in the 
model realized in CalcHEP/CompHEP packages which we describe below. 

The interactions with gauge bosons are described by , obey the SU{3)c x SU (2) x 
U{1)y Standard Model symmetry. 
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Table 1 

Eq, 



2.2,2.3 



Quantum numbers for the complete set of scalar and vector ]Q fields appearing in 
Sy, Sy, S3, R2, R2, and V^" , , f/f, C/f, f/s, respectively @. Also, TQ-lepton- 
quark couplings [Al(?9), \R{lq) and XLivq)] and notation for IQ names for the model realized in 
CalcHEP/CompHEP are presented. The particle-antiparticle convention is defined as: ^f=2 Iq 
and ^F=o — > 



The IQ — gluon interactions are described by the Lagrangian of the most general 
form for the scalar and vector Ki interactions, and Cy, respectively: 



4= E 

scalar 



(2.4) 



vectors 



Ws 



(1 - KGW^tf.HGr + ^^;ii?i^j^^a 



(2.5) 



Here, gs denotes the strong coupling constant, ta are the generators of SU{3)c, Ms{My) 
are the scalar (vector) leptoquark masses, while and Xq are the anomalous couplings 
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related to the anomalous magnetic and quadrupole moments of vector [18|. Fields ^ 
and represent scalar and vector leptoquarks, respectively. The field strength tensors of 
the gluon and vector leptoquark fields are: 



with the covariant derivative given by 



(2.6) 



(2.7) 



We omit here the analogous U(1)y x SU{2) piece of gauge interactions since it is not 
relevant for our study of IQ production at LHC, where (gauge boson - IQ) interactions are 
eventually gluon dominant. 

2.2 color factorization and model implementation into CalcHEP/CompHEP 

The appearance of vertex with 4 color particles can not be straightforwardly implemented 
into CalcHEP/CompHEP packages. The idea is to split 4-color interactions into 3-color 
vertices via the introduction of auxiliary ghosts fields. 

It is well known that for Gqq and GGG QCD vertices their color structure is factoris- 
able. For calculation of Feynman diagrams with such vertices, where color indexes can be 
convolved separately, an elegant technique is presented in |19|. However, for 4G vertex 
such a factorization is absent. To have color factorization in the case of AG vertex one can 



split this vertex into 3G vertex by means of the auxiliary tensor field Gt^,y|2C]. This field 
should have the point-like propagator 

1 



and interact with gluons according to 



(27r)4i 



SctGG 

Graphically, it can be represented as 



^ / f^,Gt^^''{x)Gf,ix)G2{x)d' 



(2. 



(2.9) 



Gi 



^<?4 G2' 



Gt 



Gs Gi 

* + Gt 

G4 G2 • 



-G3 Gi-^ 
+ 

-G4 G2^'^ 



\Gt 



^^4 



Figure 1: Splitting of a four-gluon vertex into three- gluon vertex. 

This trick was successfully implemented in CompHEP/CalcHEP packages for QCD 
Lagrangian. However in the case of the interaction of vector leptoquarks with gluon |18| 
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one needs to develop this technique further and find a general prescription for splitting 
of the arbitrary color structure. Such a prescription is described below. We consider 
separately cases of real (e.g. gluon) and complex (e.g. leptoquark) vector fields. 

In the case of a real vector field interacting with a gluon field G^, the Lagrangian 
can be presented as 

L = - du^^ + T^,|2 + /nt(T($, G), ...) (2.10) 

where T^j^y = g[^i_i, Gp], Int represents terms of gluon interaction with color matter realized 
via general derivative or the strength tensor F^^. The color index of the gluon fields is 
omitted to simplify the expressions. 

Now, we replace T^^ by an auxiliary field T^^ and confine of T = T by means of a 
Lagrange multiplier t'^'^ . We apply this procedure for the Int and T(<I>, G)^ terms, keeping 
the linear T term as it is: 

-\t^,T^^- + \t^,{T^^ - T(cl>, GY') 

+/nt(T ^ X, $,...) (2.11) 
Then, we perform the following substitution of variables 

^T^,u = {T^,^-t^^)/V2 and = i • V/^/2 (2.12) 

in order to obtain standard normalized quadratic forms for the auxiliary fields: 
L = -\\d^^u - d,^^\^ - ]^{d^^, - a,G^)T(cI>, GY" 
-h^^^^^T,^- - \^t^,^r + -J=<i>tT($,G)'^- 
+/nt(T^^/2($T-^-^t),G,...) (2.13) 

In the case of QCD interaction with scalar and fermion fields, the /nt() term does not 
contain T. Thus, the auxiliary field T is free and can be omitted. Other terms produce 
the 3 — gluon interaction, the point-like propagator (^), and the interaction of just one 
auxiliary tensor field t with gluons (^^). 

Lagrangian with a complex vector field $ also can be represented in a color-factored 
form using the same trick. In a general case, it is given as: 

L = -^\d^^u - d,'^^ + T(cl>, G)^,p + /nt(T(c^, G), T(cl>, G)*, <I>, ...) (2.14) 

Following the procedure adopted with real vector field, we introduce two auxiliary fields 
Tf^i, and tf^i, which are now complex. In terms of these fields, the Lagrangian becomes: 

L = -l\d^^. - d.^^f - ^{d^K - d,%)T{^,GY'' - ^T(ci>,G)*,(a^y^ - d'^vn 

_hr*^T,u ^ _ T(ci>, G)^^) + i(T;, - T(cl>, G)l,)t^^ 

+Int{T ^T,T{'^,Gy ^T*...) (2.15) 
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To obtain a standard normalized quadratic form of auxiliary fields, one should perform 
the following variable substitutions: 



= i ■ > butl '^t*=i- tljV2 ; (2.16) 



Note that this transformation is realized through a non- analytic manner which is legal 
for functional integrals. Now, the Lagrangian in color-factorized form can be given as 



+Int{T V2{^T - i • ^t), T($, G)* V2{^^ - i • (2.17) 
In particular, for the case of Lagrangian ( ^.51 ) one has 



T($,G)^, = -i5(G^$,-G,$^) (2.18) 

T{'^,G);, = igiG,<l>i-G,<i>l) (2.19) 
Int{T V2{'^T - i ■ ^t), T($, G)* ^/2($^ - i ■ $*)•••) = 
-ig{l - KG)^^e^^.[d^Gl - d,Gl + V2{^t1u - i^Tt^v)] 



[a^^>^ - d^"^^ + \/2($Ti; - i^t^u)] X WG" - d^G" + V2(Gt''" - iGD] (2.20) 

In the case of a complex scalar field one should introduce two auxiliary vector fields - 
and (analogous to T^^, and t^u tensors). The Lagrangian with the complex scalar 
field $ can be represented as: 

L = \d^^ + T($, G)^|2 + /nt(T($, G), T($, G)*, ...) (2.21) 

In terms of these fields, in exact analogy with the vector leptoquark field, we have 

L = + 9^$*T($, GY + T($, G)l^^'^ 

+V;V^ - vl[V^^ - T(cD, GY] - [V; - T(cl., G)l]v^^ 

+Int{T V, T($, GY V*...) (2.22) 

with the same variable substitutions: 

^v^. = {v^ - v^)/^ = {v; - vDi^ 

^^^ = i-v^,/V2 but] '^^* = i-v;jV2 (2.23) 

For Lagrangian (^^) one has: 

Ti'^,G)^ = -igG^<P (2.24) 

Ti<^,GY^ = igG^^ (2.25) 

Int{T ^ V2{^v - i ■ «>^), T($, GY ^ \/2($^ - i ■ K)...) = 0. (2.26) 
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Table 2: An example of the implementation of interactions with gluons in the 
CalcHEP/CompHEP packages for the case of vector IQ (VM) with F = -2, T3 = -1/2, Q = 
+1/3. See text for details. 



One can see that the leptoquark Lagrangians given by Eq.(^.4|)- (|2.5| ) can be rewritten 
in terms of products of only three fields. The approach described above allows us to 
represent the strength tensor as a sum of vector fields and auxiliary tensor fields. Thus, 



interactions can be presented as a trilinear fields interaction only and, therefore, have 
simple unambiguously defined color factor that can be factorized. 

As it was demonstrated above, the implementation of the complete IQ model requires 
two auxiliary tensor fields ^.t and ^.T for each IQ type. In CompHEP pO], only one tensor 



field is automatically generated, while CalcHEP [21| generates both of them. To im- 
plement the complete IQ model into CompHEP one can introduce the second tensor field 
by doubling the vector IQ fields (by introducing, say, ^> IQ fields), and use the tensor field 
of those particles {^-t) as a ^>.T field. The complete IQ model given by Eq. ( |2.2| )-( |2l5| ) has 
been implemented and tested for both CompHEP and CalcHEP packages. The complete 
models of IQ interactions relevant to LHC or Tevatron physics (i.e. models with IQ — I — q 
and IQ — gluon interactions) can be found at 

http : / /hep . pa . msu . edu/belyaev/public/pro j ects/lq/models/lq_calc . zip 
and at 

http : / /hep . pa . msu . edu/belyaev/public/pro j ects/lq/niodels/lq_comp . zip 
for CalcHEP and CompHEP, respectively. Results for IQ pair production has been com- 
pared with those of paper We agree with the results of except with the sign in 
front of \g in the expression for vector IQ pair production (for the sign convention defined 
by Eq. |2.5| ), which we found to be opposite in our calculations. 

In Table ^, we present an explicit example of the implementation of IQ interactions with 
gluons in the CalcHEP package for the case of vector IQ {VM) with F = —2, T3 = 
-1/2, Q = -M/3 (see Table |l|). Table | is a piece of the CalcHEP model of particle 
interactions. The first three fields 'PI', 'P2' and 'P3' include the names of the interacting 
particles. The last two fields 'Factor' and 'Lorentz Part' define a vertex itself. Each line of 
the table represents a particles interaction vertex. Symbols 'up', 'up.t','up.T' correspond 
to a vector IQ, first tensor IQ and second tensor IQ fields, respectively. Names with capital 
letters: 'UP', 'UP.t','UP.T' correspond to the hermitian conjugated fields. Gluon fields and 
its two tensor fields are denoted by 'G','G.t','G.T', respectively. The 'Factor' field contains 
'GG', 'MUP', 'LG' and 'Sqrt2' which denote Qs, Mup {IQ mass), Ac and \/2, respectively. 
Symbols ^mi.mj\ ^Mi.mj\ 'mi.Mj^ or 'Mi.Mj^ = 1,...3) in 'Lorentz Part' stand for 
gniimj ^ givUmj ^ gin^Mj gMiMj j^^g^j-jg tcnsors, respectively, with the indices folded with 
the respective index of vector or tensor particle in the column 'i' and 'j'. The first index 
of the tensor particle is denoted by a lower case letter, e.g., 'ml', while the second index - 
by a capital letter, e.g. 'Ml'. The index of the vector particle is always a lower case letter, 
e.g. 'ml'. The symbols 'pi.mj' or ^pi.Mf denote 'p™^ ' and 'Pj i.e. the momenta of the 
interacting particles with the respective Lorentz indices. The symbols 'pi.pj' denote the 
momenta product, Pi-pj, of two particles. The symbol 'KG' in the 'Lorentz Part' denotes 
the anomalous coupling kg- A detailed explanation on the implementation of models of 
particle interactions can be found in the CompHEP manual [pC|]. 

3. Signal rates 

Leptoquarks can be produced in quark-quark, quark-gluon and gluon-gluon interactions. 
Their production at the LHC is dominated by the strong interaction of IQ with gluons. 
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Figure 2: Diagrams for pair production in quark-quark and gluon-gluon fusion. 
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Figure 3: Diagrams for IQ single production in quark-gluon fusion 

and one can safely neglect IQ interactions with photons, Z and bosons. In the case of 
quark-quark and gluon-gluon fusion, IQ produced in pairs are shown by Feynman diagrams 
in Fig. ||. The IQ pair production cross section is defined by IQ mass and coupling of the 
strong interactions - Ofg. In the case of vector IQ pair production the total cross section 
also depends on the anomalous and Ac couplings. 

In the case of quark-gluon fusion, IQ are produced singly in association with leptons. 
In comparison to IQ pair production, the production rate of IQ single rate depends also on 
IQ — quark — lepton Yukawa type coupling appearing in Eq. (^)-( [2l^ ). The corresponding 
diagrams are shown in Fig. |3|. 

Both processes — IQ single and pair production — give the same striking signatures: 

1. 21 + jets when IQs decay into a lepton and a quark; 

2. / + jets+ for both processes when singly produced IQ decays into a neutrino and 
a quark and one of the IQ produced in pair production, decays into a lepton and a quark 
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while the other decays into a neutrino and a quark. The same signature will take place 
for IQ single production if IQ is produced in association with a neutrino but decays into a 
lepton and a quark. 

3. jets+ signature when IQs decay into a neutrino and a quark and single IQ is 
produced in association with a neutrino. 

One of the important message which we would like to convey in this article is that 
single and pair productions must be considered together. Because of similarity in 
their signatures, it is impossible to separate signals from IQ single and pair production 
completely. 

In our study, for numerical calculations we have used the CTEQ6L set for parton 
distribution function (PDF) [p4| and chosen the QCD scale equal to the IQ mass. There 
are many scalar and vector IQ species with different isospin, charge and IQ — quark — 
lepton couplings (as given in Table ||), however, the total IQ decay width depends on one 
parameter (given that IQ mass and as are fixed), the sum of IQ — quark — lepton couplings 
squared, A^^^ = \\{lq) + \\{lq) + X\{vq') (see Table |l| and Eq. (^. Without loosing 



generality, for our analysis, we have chosen scalar(5'0) and vector (VM) IQ (see Table |l]), 
which couples to both - the charged lepton and a quark as well as to a neutrino and a 
quark. For the reference point, we have chosen Ag// to be equal to the electromagnetic 
coupling, e = V'ivra ~ 0.312 for a = 1/128. Choosing Ag// completely defines IQ decay 
width: 

50with F = -2,r3= 0, Q = +l/3has r = M^^^iM±^iM), (3.1) 

()<? (la) + A^ (va') 

while with F = -2, Tg = -1/2, Q = +1/3 has F = Mjn ^ \ 'i J _ ^32) 
Both, 5*0 and VM, decay into ne^ and di> with 50% branching fraction for each decay 



channel for our choice of parameters. For Mjq = 1 TeV, y^l^. + ^SL ~ \J 92R~^ dii ~ 
Xeff = e, the total width for SO is r(S'O) = 1.95 GeV while r{VM) = 1.30 GeV. 

In Fig. ^ we present the total cross section of scalar IQ production as a function of 
IQ mass for IQ pair (left) and single (right) production. One should note, that single IQ 
and single IQ production rates at the LHC are drastically different since the first of them 
is initiated by valence quarks while the other one ~ by sea quarks. Below, we quote the 
sum of the cross sections for IQ and IQ single production. IQ single production rates also 
depend on to which quark(s) it couples. In Fig. |^ one can see the difference in rates for 
single IQ production when it is produced in association with charged lepton (higher rates) 
and when it is produced in association with a neutrino (lower rates). The apparent origin 
of this difference lays in the difference between up— and down— quark parton densities, 
respectively. We do not lose generality when studying just one type of scalar IQ, SO here, 
since the generic cross section of the scalar IQ single production could be expressed as a 
superposition of two cross sections represented by the red and black curves in Fig. Q which 
scales quadratically with the IQ — I — q coupling. The same is also true for the case of 
vector IQ production which we describe below. 
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Figure 4: Total cross section of scalar IQ (Sq) production as a function of the mass for 
pair (left) and single (right) production at the LHC. The cross section of single production 
depends on the value IQ — i ~- q coupling chosen to be y/ gi^^ + g'^^ = Ae// = e. The single 
production in association with charged lepton is represented by black curves while its production 
in association with neutrino is shown by red curves. 



For low masses (~ 100 GeV) Ki pair production process dominates over the Ki 
single production process by about one order of magnitude: ~ 1000 pb compared with 
~ 100 pb. But for larger masses (~ 1000 GeV), targeted by LHC (unless Tevatron 
will discover light IQ states earlier) the cross sections of IQ single and pair production 
become comparable. This happens because of a stronger phase space suppression of IQ 
pair production and a related faster drop of the parton (especially gluon) density functions. 
For MjQ ~ 1000 GeV, the cross section for both, single and IQ pair production, is roughly of 
the order of 10 fb. For even higher scalar IQ masses, single production starts to dominate 
over IQ pair production - for Mjq ~ 2000 GeV, the IQ single production cross section 
(~ 0.1 fb) is already one order of magnitude larger than the cross section of pair IQ 
production (~ 0.01 fb)! Of course, the cross section of IQ single production is proportional 
to the square of the IQ — e — q unknown coupling (chosen to be equal to the electromagnetic 
coupling) and can be rescaled accordingly. One can see that if this coupling is of the order 
of the electromagnetic coupling (like in the case of recent experimental limits obtained at 
HERA Q), then the scalar IQ single production should be definitely taken into account 
and combined with the studies of the IQ pair production. 

Vector IQ cross sections are presented in Fig. ^ and Fig. ^ for the cases of IQ single 
and pair production, respectively. We present results for four 'traditional' choices of kg 
and \g'- 

1. Kg = Ac = 0, Yang-Mills type coupling (YM) case 
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Figure 5: Total cross section for pair vector IQ production as a function of the IQ mass. Four 
choices of (k, A) are presented: 1) kg = = - Yang-Mills type coupling (YM); 2) k,q = 1, Xq = Q 
- Minimal coupling (MC); 3) kg = — 1, Ag = — 1 - (MM) and 4) the case of absolute minimal cross 
section (AM) in which the cross-section is minimized with respect to kc^g parameters for each 
value of M]Q. 



2. KG = 1, Ag = 0, Minimal coupling (MC) case 

3. KG = -1, Ag = -1, MM case 

4. the case of absolute minimxil cross section (AM) in which the cross section is minimized 
with respect to kg, Ag parameters for each value of Miq. 

For AM case, one can establish the absolute conservative limit on the cross section of vector 
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Figure 6: Total cross section of single vector Ej{VM) production as a function of the mass. 
The cross section of IQ single production depends on its quark content and the value IQ — i — q 
couphng chosen to be \/ sl/?, + 92L — ^eff — e. single production in association with the charged 
lepton is denoted by black curves, while its production in association with neutrino is denoted by 
red curves. Four choices of (k, A) are presented: 1) kg = = - Yang-Mills type coupling (YM); 
2) KG = 1, Ag = ~ Minimal couphng (MC); 3) kg = -1, Ag = -1 - (MM) case and 4) the case of 
absolute minimal cross section (AM) in which the cross section is minimized with respect to kg, Ag 
parameters for each value of Mjq. 

IQ since for each value of Mjq the minimal cross section is defined. One can see again, that 
the cross section for IQ single production catches up the cross section of pair production at 
^IQ — 1 TeV and starts to dominate at larger IQ masses for the same reason as in scalar IQ 
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production case. One can also notice that the cross section for MC and AM cases are quite 
close to each other for Ki masses < 1 TeV, e.g. at Mjq = 1.0 TeV they are 37.6 fb (MC) 
and 24.8f b (AM) for pair vector IQ production and 31.4 fb (MC) and 30.2 fb (AM) for 
vector IQ single production. Here, we present the sum of the cross sections, pp BQi and 
pp — > SQi', of JQ single production processes. However, at higher masses (1.5-2 TeV) 
the cross section for pair production in the AM case could be a factor of 2-3 smaller 
than for MC case. This happens because the contribution from gg IQIQ to the total 
production cross section (which is similar for MC and AM cases) vanishes for heavy IQ 
production, while the contribution from qq IQIQ starts to dominate (it is several times 
larger for MC, as compared to AM case) in this mass region. 

The cross section for YM case is typically a factor of 3-5 higher than the one for MC 
case. For example, for Mjq = 1.0 TeV one has 208 fb and 87.3 fb for IQ single and 
pair production cross section for YM case. For MM case, the cross section for pair IQs 
is 1-3 orders of magnitude higher than the one for MC case: (t^^^{MM) = 1760 fb for 
Mjq = 1.0 TeV. The MM cross section for IQ single production is about one order of 
magnitude higher, as compared to MC case: : a'^^^^'^{MM) = 284 fb for Mjq = 1.0 TeV. 

In Table ^, we summarize cross sections for IQ single and pair production for both, 
vector iVM) and scalar (50) leptoquarks at the LHC. As we mentioned above, the total 
cross section for any IQ type can be obtained from Table ^ by superimposing respective 
numbers and rescaling them for a chosen value of IQ — I — q coupling. 

In our study we use the MC set of cross sections (which is not very different from AM 
case) to establish a conservative limit on the Mjq and Ag// parameters of the model. 



4. Simulations and signal versus background analysis 

The simulations of leptoquark signal events were performed with CompHEP [^0|, the 
CompHEP-PYTHIA interface [^] and PYTHIA6.2 ||2^ program chain. The cross sec- 
tion values presented in this article were calculated using CTEQ6L parton distribution 
function (PDF) p^ . PYTHIA was used to account for initial and final state radiation and 
to perform hadronization and decay of resonances, when it was relevant. 

Besides the difference in cross section for vector and scalar IQ production one could 
expect a difference in angular correlations of vector IQ decay products compared to scalar 
IQ if the vector IQ is being produced with some polarization. We have found that, indeed, 
vector IQ is produced with non-trivial polarization at the LHC which is reflected, for 
instance, in a difference in the angular distribution between the leading pT electron and jet 
shown as an example in Fig. ^ for IQ single production case. The study of such angular 
correlation effects is, however, beyond the scope of the present article. The distributions 
of kinematical characteristics we have chosen in this study are similar for scalar and vector 
IQ production. 



The ATLFAST|25| code has been used to take into account the experimental conditions 



prevailing at LHC for the ATLAS detector. The detector concept and its physics potential 



have been presented in the Technical Proposal[26| and the Technical Design Report 



The ATLFAST code for fast detector simulations accounts for most of the detector features: 
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Table 3: Cross sections (in fb) for BQ single and pair production for vector (VM) and scalar (50) 
leptoquarks at the LHC. The cross section of single production depends on its quark content 
and the value oi IQ — £ ~ q coupling chosen to be \/ g\j^ + = ^e.ff = e. single TQ production 
in association with a neutrino is presented separately from JQ production in association with a 
charged lepton. Four choices for («;, A) are presented for vector IQ production: 1) acq = Ag = 
- Yang-Mills type coupling (YM); 2) kg = 1, Ac = - Minimal coupling (MC); 3) kg = -1, 
Ag = — 1 - (MM) case and 4) the case of absolute minimal cross section (AM) in which the cross 
section is minimized with respect to kg, Ag parameters for each value of Mjq. 



jet reconstruction in the calorimeters, momentum/energy smearing for leptons and photons, 
magnetic field effects and missing transverse energy. It provides a list of reconstructed jets, 
isolated leptons and photons. In most cases, the detector-dependent parameters were 
tuned to values expected for the performance of the ATLAS detector from full simulation. 
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The electromagnetic calorime- 
ters were used to reconstruct the 
energy of leptons in cells of dimen- 
sions ArixA(f) = 0.025x0.025 within 
the pseudorapidity (77) range —2.5 < 
rj < 2.5; (p is the azimuthal angle. 
The electromagnetic energy reso- 
lution is given by 0.1/^/E{GeV) 0.007 
over this pseudorapidity region. The 
electromagnetic showers are identi- 
fied as leptons when they lie within 
a cone of radius AR = \J (ArfY x (A^)^ 
0.2 and possess a transverse energy 
Et > 5 GeV. Lepton isolation cri- 
teria were applied, requiring a dis- 
tance AR > 0.4 from other clus- 
ters and maximum transverse en- 
ergy deposition, Et < 10 GeV, in 
cells in a cone of radius AR = 0.2 
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Figure 7: Distributions of the cosine of the open angle 
between electron and jet in events with a scalar (solid line) 
around the direction of electron emis- and vector (dashed line) single production, 
sion. 

Jet energies were reconstructed by clustering hadronic calorimeters cells of dimensions 
Ari X A(f) = 0.1 X 0.1 within the pseudorapidity range —2.5 < r] < 2.5. The hadronic 
energy resolution of the ATLAS detector is parametrized as 0.5/a/ E{GeV)(B0.03 over this 
ij region. Hadronic showers are regarded as jets if the deposited transverse energy Et is 
greater than 15 GeV within a cone of radius AR = 0.4. 

It must be mentioned that standard parametrization in ATLFAST has been used for 
the electron resolution but detailed studies are needed, using test beam data and GEANT 
full simulation to validate the extrapolation of the resolution function to electron energies 
in the TeV range. 

In this paper, two types of signal event have been studied. 



1. Type 1(2/ + jets) signature, for which IQ produced singly in association with an 
electron, decays to an electron and a quark, while each produced in pair is required 
to decay to an electron and a quark. 



2. Type 2{l+jets + Ej!''^^^) signature, for which ]Q produced singly in association with 
a neutrino, decays to an electron and a quark. ^ For pair production, one is 
required to decay into electron and quark, while the other one - to a neutrino and a 
quark. 

^The same signature will be for single IQ production in association with the electron if ]Q decays into 
neutrino and quark. 
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Backgrounds of Type 1 signal signature are Z + jets events, where Z decays into two 
electrons, and it events where both W from the top quark decay into an electron (positron) 
and a neutrino. For the signal events of Type 2, backgrounds are VF+jets events, where W 
decays into an electron and a neutrino, and it events where one W decays into an electron 
and neutrino, and the second W decays into jets. 

Process tt(Type 1) it (Type 2) Z + jet W + jets 
a (pb) 6.6 11 665 17. 

Table 4: Cross section (7(pb) from PYTHIA for the backgrounds after pre-selection cuts. 

In order to enrich the event statistics in the region of high invariant masses, simulated 
background events were pre-selected in PYTHIA for hard 2 — > 2 process with transverse 
momentum, pT > 200 GeV (100 GeV for Z + jet events), defined in the rest frame of 
the hard interaction, followed by the standard initial and final state radiation technique. 
Top quark pair production and W + jets background events for Type 2 signal events were 
additionally pre-selected with at least one electron and > 100 GeV. Corresponding cross 
sections for the backgrounds, used for this article, are shown in Table ^ for the events, 
passing pre-selection cuts. Since the signal events would consist of events originated from 
the single or pair production of leptoquarks, we tried to define a unique set of cuts effective 
for the background suppression against combined IQ pair-l-single signal. 

4.1 Type 1 signal events 

Here, we present the analysis of IQ events of Type 1. The signal signatures consist of at 
least one jet and two electrons. 

Due to their large mass, leptoquarks would produce events with large transverse energy. 
The scalar sum of transverse momenta of all charged particles in the event, Ht, is presented 
for single, Fig. ^(left), and pair. Fig. ^(right), scalar production. As can be seen, an 
appropriate choice of a Ht value should suppress the bulk of background events. 

The following cuts were used to separate signal from background: 

• The transverse momentum of two electrons were required to be at least 90 GeV (100 
GeV in the case of M^j > 750 GeV). 

• At least one jet was required with minimum transverse momenta of 70 GeV (90 GeV 
in the case of Mjq = 750 GeV and > 100 GeV/c for higher masses). 

• The invariant mass of two electrons was required to be larger than 150 GeV in order 
to veto the dominant Z + jets background. 

• Events with at least one b-jet were vetoed to suppress ti background. 

• The scalar transverse momentum sum of all selected particles in the event, Ht, was 
required to be at least 800 GeV. 
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Figure 8: Distribution of the Ht variable (see text) for events of single (left) and pair (right) 
production of scalar leptoquarks and corresponding backgrounds. All distributions are normalized 
to unity. 

The resulting invariant mass distributions for the electron-jet system are presented in 
Figs. 1^ for combined single and pair production of leptoquarks of various masses. Signal 
distributions are presented for scalar (hatched area) and vector leptoquarks. The dashed 
line shows the pair production contribution to the total signal spectrum. The contribution 
of the single production to the total spectrum gradually increases with the IQ mass. The 
cross section for vector leptoquarks corresponds to MC case {kq = 1, Ac = 0). 

All possible mass combinations between two leading pT jets and electrons are allowed, 
which leads to some broadening of the signal distributions for the single production case. 

As can be seen from these figures, the Z + jets background is essentially dominant. 
The pair production of top quarks gives a rather small contribution and can be seen only 
for the lowest leptoquark mass case studied. 

The signal statistical significances are given in Table H. The data are presented for 
combined (single-|--fi5 P^'i^' production) signal efficiencies, with the number of events shown 
for the total background and for IQ single and pair production for different masses of 
scalar and vector leptoquarks. The data are presented for an integrated luminosity of 
L = 300/6"^ 

4.2 Type 2 signal events 

In this section we present the analysis of signal events with an electron, jets and a neutrino. 
The signal signature is at least one jet, electron and missing transverse momenta. The 
scalar sum of transverse momenta of all charged particles in the event, Ht, is presented for 
single. Fig. |l^ (left), and pair production. Fig. |l^ (right) of scalar K^. As can be seen, an 
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Figure 9: Invariant mass distribution of electron and jet for events of single and pair production 
of scalar and vector leptoquarks mass 500, 750, 1000 and 1250 GeV for Type 1 signal events. 



appropriate choice of a Ht value again should suppress the bulk of background events. The 
presence of a neutrino in the IQ signal events of this type, provides a possibility to suppress 
relevant backgrounds using a missing transverse momentum variable. In Fig. [l^ (right), 
we plot the orbital (0)-angle difference, At/), between electron and the missing transverse 
momentum vector. Distributions are shown for scalar and vector leptoquark decays and 
corresponding backgrounds. The background spectra reveal the emission of an electron and 
a neutrino along the same direction, which is reasonable, since in background events those 
particles are produced in W decays. The signal distributions bear mostly the back-to-back 
emission feature. The following set of cuts was worked out to effectively separate the 



-19- 



IQ Mass 

y\j*a V J 


Background 


Scalar 


Vector 




Tig-Pair 


Stot/VB 


IQ-Smg\e 


Tig-Pair 


Static 


500 


614 


1835 


14529 


660 


25907 


155374 


7316 


750 


264 


335 


1033 


84 


2902 


8571 


706 


1000 


132 


54 


112 


14 


392 


894 


112 


1250 


40 


12 


17 


5 


70 


140 


33 


1500 


28 


4 


3 


1 


16 


20 


7 



Table 5: Number of signal events for IQ single and pair production versus total background and 
respective significance Stot I VB for combined signal {K^ pair -I- single production) . Results are for 
an integrated luminosity oi L = 300/6^^ and Type 1 signal signature. 




Hj, GeV Hp GeV 



Figure 10: Distribution of Ht variable (see text) for events of single (left) and pair (right) produc- 
tion of leptoquarks and corresponding backgrounds. Decays of include electrons and neutrino. 
All distributions are normalized to unity. 



signal from background: 

• The transverse momentum of the electron was required to be at least 100 GeV. 

• At least one jet was required with a minimum transverse momentum of at least 100 
GeV. 

• The transverse mass of an electron and the missing transverse momentum vector 
(P™**, see Fig. 12) was required to be larger than 200 GeV in order to veto the 
dominant W + jets background. 

• Events with at least one b-jet were vetoed to suppress ti background. 
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• The scalar sum of transverse momenta of all selected particles in the event, Hx, was 
required to be at least 400,(600, 800, 1000, 1200) GeV for IQ masses of 500, (750, 
1000, 1250, 1500) GeV, respectively. 

• The orbital angle difference between an electron and the transverse missing momen- 
tum vector was required to be greater than 0.8 radians. 
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Figure 11: The distribution of the orbital (/)-angle difference between an electron and the missing 
transverse momentum vector, A0, in radians, for signal (scalar and vector IQ) and background 
W + jets and ti events, single (left) and pair (right) IQ production are presented. 

While the angular spectra of scalar and vector IQ of 500 GeV mass are quite close to the 
background spectra in the forward hemisphere, the signal distributions for larger masses 
allow the separation of the signal from the background. The resulting invariant mass 
distributions for the electron-jet system of the signal events of Type 2 are presented in 
Fig. |l3|for the IQ mass of 500 GeV (left side) and 750 GeV (right side). Signal distributions 
are presented for scalar leptoquarks (hatched area) and vector leptoquarks for the minimal 
coupling set. The dashed line shows the pair production contribution to the total signal 
spectrum. Similarly to the Type 1 signal events, the contribution of the single production to 
the total spectrum gradually increases with the IQ mass. The signal statistical significances 
are reported in Table |^ combined for IQ single and pair production for different IQ masses. 
Table ^ also presents number of events, separately for for IQ single and pair production 
as well as for background events. The data correspond to an integrated luminosity of 
L = 300/6-^ 



4.3 IQ mass reach of LHC 



The IQ mass reach of LHC is shown in Fig. 14 for combined single and pair production 
processes and two types of signal signatures. The results for scalar and vector IQ are 
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Figure 12: The distribution of the transverse mass of an electron and the P™*^* vector for signal 
(scalar and vector TQ) and W +jets and tt backgrounds, single (left) and pair (right) production 
are presented. 
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Table 6: Number of signal events for single and IQ pair production versus total background and 
respective significance Stot/VB for combined signal (pair + IQ single production). Results are for 
an integrated luminosity oi L = 300/6^^ and Type 2 signature. 



presented separately. The data are presented for an integrated luminosity of L = 300/6~ . 
One can see that scalar leptoquarks can be accessible at LHC up to masses < 1.2 TeV 
while vector IQ for MC case can be discovered for masses < 1.5 TeV. LHC can exclude 
IQ at 95% CL with masses about 200 GeV above 5a discovery IQ limit, i.e. with masses 
< 1.4 TeV and < 1.7 TeV for scalar and vector (MC) IQ, respectively. 

Let us remind, that we have chosen IQ — I — q coupling Ag// = e and contribution 
from IQ single production rescales quadratically with this coupling. For other values of 
Ae//, the new LHC reach can be easily found by using LHC reach Tables 5 and 6. One 
can see that Type 1 signature {2i + jets) looks more promising compared to Type 2 events 
{2£ + jets+:pT)- 

One should notice that, for the chosen Ag//, the contribution from IQ single production 
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Figure 13: Invariant mass distribution of electron and jet for events of single and pair production 
of scalar and vector leptoquarks for m — 500 GeV (upper left) and m — 750 GeV (upper right) and 
TO = 1000 GeV (below) for Type 2 signature. 



is 30-50% to the total number of signal events for the IQ mass at the discovery limit. 
Therefore, single and pair production should be studied together at the LHC ^. The 
complementarity of single and pair IQ channels is clearly illustrated in our final Fig. 14, 
where we present scalar and vector leptoquark mass reach of the LHC in the {Mjq — \eff) 

^Eventually, the case of single production of of the second and third generations is qualitatively 
different. In this case only pair IQ would give the major contribution to the signal rates, unless \eff is 
too large (which might be not allowed by other experimental constraints) to enhance IQ single production, 
suppressed due to initial sea-quarks PDFs. 
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Figure 14: Leptoquark mass reach of the LHC at 5a level for Type 1 and Type 2 signal signatures 
for scalar (solid line) and vector (dashed line) for an integrated luminosity of L = 300/5~^. For 
vector SQ, the MC choice has been selected. The contributions of single and pair signal are 
combined. 

plane. For example, for Ag// = 1, IQ single production allows the extension of LHC reach 
for MiQ by about 400 GeV! 

5. Conclusions 

In this paper, we present a new detailed study of IQ production and decay at LHC at the 
level of detector simulation. 

We treated IQ single and pair production together and have worked out a set of 
kinematical cuts to maximize significance for combined single and pair production events. 

It was shown, that combination of signatures from IQ single and pair production not 
only significantly increases the LHC reach, but also allows us to give the correct signal 
interpretation. Our results are summarized in Tables 5 and 6 and Figs. 14 and 15. In 
particular, the LHC can discover IQ with a mass up to 1.2 TeV and 1.5 TeV for the case 
of scalar and vector IQ, respectively (for Ag// = e), and IQ single production contributes 
30-50% to the total signal rate. 

In this work, the most general form of scalar and vector IQ interactions with quarks 
and gluons has been implemented into CalcHEP/CompHEP packages, which was one of 
the primary aspects of this study. 
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Mass, TeV 

Figure 15: Scalar and vector leptoquark mass reach of the LHC for {Miq A) plane for pah- 
production (dashed line) and for combined single and pair production (solid line) . Type 1 signal 
reach contour is denoted by stars while Type 2 signal reach contour is marked by solid dots. 
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